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The phase parameters of matched-fihering searches for continuous gravitational-wave signals are 
sky position, frequency and frequency time-derivatives. The space of these parameters features 
strong global correlations in the optimal detection statistic. For observation times smaller than one 
year, the orbital motion of the Earth leads to a family of global-correlation equations which describes 
the "global maximum structure" of the detection statistic. The solution to each of these equations 
is a different hypersurface in parameter space. The expected detection statistic is maximal at the 
intersection of these hyp ersurf aces. The global maximum structure of the detection statistic from 
stationary instrumental-noise artifacts is also described by the global-correlation equations. This 
permits the construction of a veto method which excludes false candidate events. 



PACS numbers: 04.80.Nn, 95.55.Ym, 95.75.-z, 97.60. Gb 



I. INTRODUCTION 

The emission of continuous gravitational waves (CW) 
is expected, for instance, from spinning neutron stars 
with non-axisymmetric deformations. If the system is 
isolated, it is losing angular momentum through radia- 
tion and is slowing down. Therefore the gravitational- 
wave frequency would be slowly decreasing for this long- 
lasting type of signal. Such CW sources are among the 
primary targets of Earth-based, laser-intcrferometric and 
resonant-bar detectors. 

The terrestrial location of the detectors generates a 
Dopplcr modulation of the signal caused by the detector's 
motion relative to the solar system barycenter (SSB). 
The observed phase therefore depends on "phase param- 
eters" , which describe the intrinsic frequency evolution 
and the source's sky location. In addition, there is a 
time-varying amplitude modulation due to the antenna 
patterns changing with the Earth's spinning motion. The 
latter variations depend on the "amplitude parameters" , 
which are the two polarization amplitudes, and the po- 
larization angle of the gravitational wave. 

To extract CW signals buried in the detector noise, 
the optimal data analysis scheme is derived in [l| based 
on the principle of maximum likelihood detection lead- 
ing to coherent matched filtering. It is shown, that the 
amplitude parameters together with the initial-phase pa- 
rameter can be eliminated by analytically maximizing the 
detection statistic, such that the search space is just the 
phase parameters: sky position, frequency and frequency 
time-derivatives. This detection statistic is commonly 
referred to as the ^-statistic. For a given sequence of 
data, wide-band all-sky searches evaluate the JF-statistic 
over a large number of template-grid points in param- 
eter space. The parameters of the templates for which 
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a predetermined threshold is exceeded are registered as 
candidate events for potential gravitational-wave signals. 

In the local parameter-space neighborhood of a given 
signal one can define a metric [3,11, 13] from the frac- 
tional loss in expected JF-statistic. This fractional loss 
defines the dimensionless "mismatch" fi. Let A define a 
vector of phase parameters for a template. If As denotes 
the signal's phase parameters, then the metric fl'ij(As) is 
obtained by Taylor-expanding the mismatch at As with 
respect to the small parameter offsets AA = As — A to 
quadratic order: = J2ij 5u(^s) AA^AA-? +0[{AXf]. 

This work identifies the global parameter-space regions 
where the detection statistic T is expected to have large 
values close to maximal without restriction to the local 
neighborhood of the signal location. In this paper, these 
regions are referred to as the ^^global large-value struc- 
ture^ . To find the global large- value structure of the !F- 
statistic, a simplified detection statistic J^* (approximat- 
ing J-') is considered. The locations in parameter space 
where J-* is expected to be maximal are referred to as 
the global maximum structure^'' . For increasing parame- 
ter offsets from the given signal's parameters this global 
maximum structure of T* (and therefore also the global 
large-value structure of T) is found to become signifi- 
cantly different from local approximation obtained from 
the metric. 

A previous study examined monochromatic signals 
in the restricted phase-parameter space of sky location 
and frequency. It was shown that in some different ap- 
proximate detection statistic, such signals can generate 
"circles in the sky" while searching a range of template 
frequencies. The collection of these circles forms a two- 
dimensional hypersurface in the three-dimensional pa- 
rameter space, described by a single equation. Here, even 
in this restricted parameter space, this hypersurface is 
only an approximation of the description of the global 
maximum structure of the .F*-statistic to first order in 
observation time T . 

The present work shows that the global maximum 
structure of the detection statistic T* is described by 




FIG. 1: Schematic drawing of intersection of two representa- 
tive global-correlation hypersurfaces Tii and Ti.2- The region 
of intersection, shown by the solid black curve, describes the 
global maximum structure of the simplified detection statis- 
tic JT*, because the expectation value of JT* is maximal in 
these locations. For illustration purposes only two hypersur- 
faces are shown, whereas in general one might need to consider 
more hypersurfaces. 



a separate equation for each order of T. The solu- 
tion to each of these equations is a different hypersur- 
face in parameter space. Therefore, it is this family 
of global- correlation hypersurfaces, which describes the 
global maximum structure of J-'*: The detection statis- 
tic J^* is expected to be maximal at the intersection of 
these hypersurfaces. This idea is illustrated schemati- 
cally in Figure [TJ The same results also apply when con- 
sidering the generalization to non-monochromatic signals 
allowing for an intrinsic frequency evolution of the source. 

This paper is organized as follows. Section |TT] briefly 
reviews the matched-filtering method for CW sources. 
Section [IIII describes an approximate signal model which 
leads to the simplified detection statistic allowing the 
analytical exploration of its global maximum structure. 
Section IIVI presents the global-correlation equations of 
the phase-parameter space and illustrates the geometry 
of the global-correlation hypersurfaces. In addition, the 
search-parameter regions are derived for which the ap- 
proximations used are valid. In Section |V] the predic- 
tions made by the global-correlation hypersurfaces are 
compared to a fully coherent JF-statistic search in data 
sets of software-simulated sources with no detector noise 
as well as of hardware-injected CW signals in the pres- 
ence of real detector noise. In Section IVII the global 
large-value of the detection statistic which is caused by 
stationary instrumental-noise artifacts mimicking astro- 
physical sources is found to be also well described by 
the global-correlation equations, and thus a veto method 
is constructed. Section IVIII discusses the effects of the 
Earth's spinning motion in the context of the present 
topic. Finally a concluding section follows. 



II. OPTIMAL FILTERING FOR CONTINUOUS 
GRAVITATIONAL- WAVE SIGNALS 

In the absence of any signal, the detector output data 
time series x(t) at detector time t only contains noise 
n{t), which is assumed to be a zero-mean, stationary and 
Gaussian random process. If a signal h(t) is present, 
the noise is assumed to be additive, x{t) = n{t) + h(t). 
Denote by tssB the time measured at the solar system 
baryccntcr. For a detector at fixed position and orienta- 
tion, at the SSB the continuous gravitational- wave signal 
is described by a sinusoid of constant amplitude and a 
phase given by 



(fc + 1)! 



j-k+l 



(1) 



where $o is the initial phase, /'■*'■' = / denotes the fre- 
quency, and j('^>o) the fc'th frequency time-derivative, 
defining every parameter j'*^' at t = at the SSB. The 
integer s > denotes the number of frequency time- 
derivatives to be taken into account, therefore it is set 
case of an isolated rapidly rotating 
neutron star with non-axisymmetric deformations and 
negligible proper motion Q, the waveforms correspond- 
ing to the plus (+) and cross (x) polarizations are 

/i+(t) = ^+ sin^'(i), /ix(t) = cos^'(<). (2) 

The Earth's motion with respect to the SSB leads to 
Doppler effects causing amplitude and phase modulations 
of the CW signal received at the detector. Define n as 
the constant unit vector pointing from the SSB to the 
source. Neglecting relativistic and higher order correc- 
tions, a wavcfront arriving at the detector at time t, 
passes the SSB at time 



t + 



r{t) ■ n 



(3) 



where the vector r{t) connects from the SSB to the de- 
tector, and c is the speed of light. 

It is shown in jlj that the resulting phase evolution 
of the continuous gravitational- wave signal can be repro- 
duced without significant loss in signal-to-noise ratio by 
the model 
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The received signal is also amplitude modulated by the 
varying antenna-pattern of the detector due to its motion 
with the rotation of the Earth. The dimensionlcss signal 
response function h{t) of an interferometric detector to 
a weak plane gravitational wave in the long-wavelength 
approximation is a linear combination of the form: 



hit)=F+{t)h+{t)+F^{t)h^it), 



(5) 
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where -F+,x are caUed the antenna-pattern functions, re- 
sulting in the amplitude modulations from Earth's spin- 
ning motion. They lie in the range — 1 < -F+,x < 1, and 
depend on the orientation of the detector and source, and 
on the polarization angle ip of the waves. 

The optimal detection statistic [H, 01 obtained from the 
likelihood ratio A defines the matched filter 
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lnA= — 



{x\\h)-^{h\\h) 



(6) 



where Sh is the one-sided noise strain spectral density 
which is assumed to be constant over the narrow band- 
width of the signal, and the inner product is defined as 



J J-T/2 



x{t) y{t) At, 



centering the coherent observation-time interval of dura- 
tion T around i = 0. Replacing the amplitude parame- 
ters in Equation ([6]) by their values which maximize In A, 
the so-called maximum likelihood (ML) estimators, de- 
fines the detection statistic T as 



= InA 



ML- 



(8) 



III. MATCHED-FILTERING DETECTION 
STATISTIC OF A SIMPLIFIED SIGNAL MODEL 

A. The simplified signal model 

The phase of the continuous gravitational-wave signal 
is expected to change very rapidly at the detector site 
on the Earth over a characteristic time length of typi- 
cally less than ten seconds, whereas the amplitude of the 
signal changes with a period of one sidereal day. As a 
result Q, the detection of a CW signal requires an ac- 
curate model of its phase, because even 1/4 of a cycle 
difference between template and signal can lead to a loss 
in signal-to-noise ratio of 10%. Whereas modeling of its 
amplitude is less critical. Therefore, the antenna pattern 
functions F+.x are assumed to be constant, so that the 
signal model ^ becomes 



h{t) = Ai cos $(t) + A2 sin $(t) , 



(9) 



where Ai_2 are defined to be the constant effective ampli- 
tudes. The validity of this approximation is investigated 
using Monte Carlo simulations in 

The vector r{t) connecting the SSB and the detector 
can be decomposed into an orbital component rorb (t) and 
a spin component rspin(t) as 



r{t) 



I'orb 
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spin 



it), 
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where rorb(0 represents the vector from the SSB to the 
Earth's barycenter, and rspin(i) is the vector from the 
Earth's barycenter to the detector. Thus, substituting 
the decomposition (fTU|) into Equation ^ one can write 



separately the orbital component (^orb(0 the spin 
component ^spin(i) in the phase model Q as 



\k=0 



$(0 = 27T \^J2^^I^t''+^ J +^orb(t)+<^spi„(i), (11) 

where 



M)^2,'-^(±!^A. ,12, 

\fc=o ■ / 

The orbital motion of the Earth has a period of 
one year, so its angular frequency is fiorb = 27r/lyr. 
Fully coherent all-sky searches for observation times T 
much larger than a few days are computationally pro- 
hibitive [l|, Thus for computationally feasible co- 
herent searches or coherent stages of hierarchical mul- 
tistage searches 0, [l^ the typical observation time base- 
line would be of order a few days. Therefore, only ob- 
servation times T are considered, which are much shorter 
than one year: ^lorbT ^ 2it. Then the component rorb(i) 
will vary slowly and one may use a Taylor expansion at 
time to with 



i=0 



it - toY 



(14) 



where rdrb(^o) denotes the ^'th derivative with respect to 
time of rorb(i) evaluated at t^. Without loss of generality 
we may choose to = as the midpoint of the observation 
of duration T in the following discussion. Define ^ = 
rorb(0)/c with ^ = 1^1 such that 



M) 



(0) 



and 



(15) 



Together with Equation the orbital component of 
the phase p2)l can be written as 
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(16) 



The spinning motion of the Earth has a period of one 
sidereal day (Isd), which translates into an angular fre- 
quency of f^spin = 27r/lsd. The corresponding aver- 
age velocity of Vspin/c ~ 10"^ is two orders of mag- 
nitude smaller than the corresponding orbital velocity, 
Voih/c ~ 10^"'. In what follows we neglect the contri- 
bution of the spin component </)spin(t) to the phase (|lip . 
Section IVIII will discuss in detail the effects of the spin 
component 0spin(O in terms of the matched- filtering am- 
plitude. 
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Using 4>oT:h(t) in the form of Equation (|16p . we refer to be written as 
^orb (t) as the "orbital phase model" : 



$orb(0 = 27r 



(18) 
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where the coefficients Um of the power series are defined 
By rcparameterization the orbital phase model (|17p can by 
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(19a) 
(19b) 

(19c) 
(19d) 



r 



so that for arbitrary order of m > the coefficient Um is 
obtained as 



f=0 



£!(to-£)! 



,(m-f) 



(19e) 



B. The simplified matched-filtering detection 
statistic 



By analogy to Equation ^ , we refer to In A* as the log 
likelihood function of the simplified CW signal model de- 
scribed in previous Section IIII Al Maximization of In A* 
with respect to the unknown amplitudes A12 yields their 
ML estimators. By substituting the ML estimators back 
into In A*, the simplified detection statistic J-* is defined 
as 

T ,, ,,2 T 



(.T||e— ) 



(20) 



2 Sh 2 Sh 

where the detection-statistic amplitude X has been de- 
fined through 

A'^(x||e-^*-), (21) 

using the orbital phase model <I>orb- 

For further simplicity consider a data set x{t) which 
only contains unit-amplitude signal s{t), such that 

x{t) =n[s{t)]. (22) 

Let the phase-parameter vector As = ({/i'^''}, ng) define 
the phase ^^^^(i) of the signal. Then s{t) can be ex- 
pressed as 



sit) 



and one obtains x{t) — cos^^^^^{t). The difference in 
phase A$orb(i) between the phase •I'orblO ^^'^ signal 
and the phase $orb(i) of a template A = ({/^''■^},n) is 
defined by 



A$orb(t) = Krhit) - $orb(i). 



(24) 



The maximization of J-* is equivalent to maximiz- 
ing I A" p. Using Equation (j23|l one may rewrite the sim- 
plified matched-filtering amplitude X as 



T 

1 
T 



T/2 



'T/2 
T/2 



T/2 



x{t) e-**"''^ dt 



.b(t)+e-4ej*)+*o,b(t)] di.(25) 



Dropping the rapidly oscillating term in Equation 
yields 



T 



T/2 



T/2 



(26) 



Thus, Equation p6|) shows that \X\ has a global maxi- 
mum of I A" I = 1; if during the observation time interval T 
the phase difference A$orb(i) is stationary: 



aA$o,b(t) 
dt 



= 0. 



(27) 



Later, we will demonstrate that for the orbital phase 
model, unity of \X\ can be obtained not only for the 
case where all individual phase parameters exactly match 



(23) (A/W = /, 



(k) 



0, An 



n = 0) and 
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so A<i>orb(i) = 0, but also for different non-zero offsets 
(A/^'"'^ ^ 0, An ^ 0) which compensate each other to 
achieve dA^orh{t)/dt w and thus lead to a value of jA"! 
close to unity. 



IV. GLOBAL-CORRELATION 
HYPERSURFACES IN PARAMETER SPACE 

A. The global-correlation equations 

The central goal of this work is to identify those loca- 
tions in parameter space where the simplified detection 
statistic J^* is maximal, which corresponds to regions 
where is one. Consistent with Equations ^8]) . 
and ((24|) . the differences between the coefHcients u,n of 
the template's phase and the coefficients of the sig- 
nal's phase are defined by 



Am,, 



Thus, Equation can be expressed as 



(28) 



(29) 



and finally one can rewrite Equation (|25p as 



X = 



e 



T 



lAtio rT/2 / oo \ 

/ exp i ^ Au,n dt. (30) 



It is apparent that \X\ does not depend upon the zero- 
order term Auq, and therefore the same holds for the 
detection statistic T*. The values of Aum for which \X\ 
attains its maximum of \X\ = 1 consistent with (|27|1 arc 
obvious from Equation (|29p . namely when 



gA$o,b(t) 
dt 



= m Au„i t^' 

m=0 



0, 



(31) 



the following central result is obtained: the family of 
global- correlation equations which describes the global 
maximum structure of J-* , 



0, 



(32) 



where m > 0, because J-* is independent of Auq. Be- 
cause (l,t,t^, . . . ) is a basis of the vector space of real 
polynomials, the zero vector can only be represented by 
the trivial linear combination, which is given by Equa- 
tion 

Thus, the first global-correlation equation given by 
Alii = can be rewritten as 



(1) 



(33) 



where Ki = u\/2'k represents a constant defined by the 
signal's phase parameters. As a side remark, it should 
be mentioned that Equation p3p is a generalization of 



the first-order global-correlation equation found in Q to 
signals with non-zero frequency time-derivatives. But 
in order to describe qualitatively the global maximum 
structure of the simplified detection statistic in parame- 
ter space, considering only the first order (as done in Q) 
might not be sufficient. As a matter of fact, this is shown 
in Section IIVDI and will be confirmed by analyzing sim- 
ulated signals later on. 

Therefore, continuing to the next order in time, one 
obtains from the second condition Au2 = the following 
relation, 

/(I) + / 1 (2) . ^ ^ 2/(1) I . n + ^ n = (34) 

where the constant K2 = u^/^-k is defined by the signal's 
phase parameters. 

In general, this scheme can extended to arbitrary order 
m > 0. Thus, the family of global-correlation equations 
represented by Aum = can be written in the form 



^ + / 



n 



i\{m-t)\ 



(35) 



where the signal's phase parameters determine the con- 
stant K,n = w^j/27r. In the following the geometry of the 
solution to the family of equations given by ()35p will be 
explored and the hypersurfaces they describe in parame- 
ter space will be illustrated. 



B. The geometry of the global-correlation 
equations 

Let the space of phase parameters A = ({/C^)},!!) be 
a manifold denoted by V. Previously we defined n as a 
unit vector pointing to the source's sky location in the 
SSB frame of reference. A position on the sky can be 
determined by two independent coordinates, for exam- 
ple one can use equatorial coordinates of right ascension 
(RA) and declination, denoted by a and S, respectively. 
In these coordinates: n = (cos 5 cos a, cos 5 sin a, sin 5). 
Recall that the integer s is related to the number of 
spin-down parameters considered (/('=>^) = 0), thus the 
phase-parameter space dimensionality is s + 3. 

By inspection of Equation (|35p . it is obvious that for a 
given continuous gravitational-wave signal the set of so- 
lutions to each global-correlation equation Aum = 0, is a 
hypersurface in the search parameter space V. Denoting 
each hypersurface by Tim one may write 



H,n = {XeV : Au„, = 0}, 



(36) 



given the signal's phase parameters u^. The dimension- 
ality of each hypersurface Tim is s -\- 2. 

Using the parameters u,„ of Equation ^9)) define the 
vector u = {{um}) = {ui,U2, ■ ■ ■) and denote the corre- 
sponding parameter-space manifold by U. Let the vec- 
tor gU denote the signal's parameters, and define the 
difference Au = u^ — u = (Ami, A?i2, . . . ). Considering 
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the detection-statistic amplitude \X\ as a function of Au, 
then this function is extremal with respect to the param- 
eter Um along the hypersurfaee Hm defined by Aum ~ 0, 



d |A'(Au) 



dUr, 



0. 



(37) 



Ati,„=0 



Obviously, on the intersection of these hypersurfaces Tim 
at Au = 0, \X\ is maximum with respect to all the 
parameters u„i: V |A'(Au = 0)| = 0. Thus, it follows 
that the simplified detection statistic J^* is also expected 
to be maximal at the intersection of these hypersur- 
faces. Therefore, the global maximum structure is de- 
fined by the intersection of the global-correlation hyper- 
surfaces Tim- This idea has been shown schematically in 
Figure [TJ In the absence of noise, all candidate events 
produced by a given CW signal which follows the simpli- 
fied model introduced in Section UlI Al will be located on 
the hypersurfaee Tii, described by Equation (j33p . From 
all candidate events located on TCi those will belong to 
loudest ones (have largest values of J-"*) which are lo- 
cated at the intersection with hypersurfaee Ti.2, which is 
described by Equation ([M)) . For each higher order m this 
behavior carries itself forward in the same way. 



C. A visualizing example of the global-correlation 
hypersurfaces 

As an illustrative example visualization of the geo- 
metrical structure formed by the global-correlation hy- 



persurfaces (j35j), we choose the four-dimensional phase- 
parameter space (s = 1). Thus, the four dimensions are 
sky position (right ascension a and declination 6), the 
frequency /, and the first spin-down parameter f^^\ For 
demonstration purposes, consider an exemplary contin- 
uous gravitational-wave signal with the phase parame- 
ters /s = 100.0 Hz, as = 2.0 rad, Ss = -1.0 rad and 
/i'^ ^ -10-1" Hz/s. 

In this illustration, assume that we have a bank of 
templates available covering the whole sky, the entire 
Doppler frequency band [11|, /s ± /s x 10~^, and for 
the spin-down value of /'^^■' = — lO-^^Hz/s. The atten- 
tive reader will notice that the signal's spin-down value 
is not covered by this template value, but mismatched 
by an order of magnitude. This choice will demonstrate 
that due to the parameter-space correlations the signal 
is still expected to be detected in such a search, pro- 
ducing a detection-statistic maximum-structure which is 
predicted by the global-correlation equations. 

The global-correlation equations only for m < 2 have 
to be considered in the present example, because contri- 
butions to I A" I in (PP]) from terms beyond second order 
are negligible, as will be shown later in Section IIVDI 
There, the neglected third-order term in Equation pO]) 
is estimated and found to be irrelevant, in the case of 
the signal investigated here, and for observation times T 
approximately less than 61 hours. 

For the case of m < 2 the integral pO|) can be obtained 
analytically as follows 



X 



T 



2Ty/A 



J{Aui t+Au2 t^) 



T/2 



U2 



exp 



i Auq — 



4Au2 



erf 



(Awi 4- Au2T)c-'''/* 

2y^A^ 



- erf 



(Ami - Am2T)c- 

2y^A^ 



'ITT /4 " 



,(38) 



where the error function erf (a;) is defined by 

2 2 

eriix) = ^ e-*'di. (39) 
Jo 

Figure[2]shows \X\ given by Equation (|38|) as a function 
of the dimensionless parameters AuiT and Au2T^. This 
figure visualizes the fact that when one moves away in pa- 
rameter space from the global-correlation hypersurfaces 
at Aum = to increasing non-zero values of Aum, the 
detection-statistic amplitude \X\ decreases rapidly from 
its maximum of one towards zero, as is shown in Figure[2l 

As mentioned earlier, the two equations from the fam- 
ily of global-correlation equations (|55|) to be considered 



for the example signal studied here are: 

/ + -n + Z^D^.n = K,, (40a) 

/«+/^(^)-n + 2/W|W.n = K2, (40b) 

where in this case the constants Ki,2 arc obtained from 
the signal's phase parameters (/s,/i^\ns) as 

Ki = fs + fs^ • ns + fP I • ns, (41a) 
K2 = /i'' + /s I • ns + 2fP I • ns . (41b) 

Figure[3]presents the visualization of the hypersurfaces 
Hi and 'H2 described by Equations (|40ap and (j40bp for 
the given CW signal. As described in Scction llVBl when 
choosing s = 1, then we have dimTii dim 7^2 ~ 3. 
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FIG. 2: (Color online) Illustration of the simplified detection- 
statistic amplitude \X\ for m < 2, given by Equation (|38p . as 
a function of the dimensionless parameters AuiT and Au2T^ ■ 



In Figure [3] the subspace {/, n} is shown at the fixed 
target spin-down of /^^-' = — lO^^^Hz/s. The hyper- 
surface described by Equation (|40ap for the fixed j'"'^-' is 
depicted three-dimensionally in Figure 3(a) The 3D plot 
of Figure 



3(b) shows the hypersurface defined by Equa- 
tion (|40b() for the same signal and same template val- 



ues. Finally, Figure 3(c) combines Figures [3(a)] and [3(b)] 
showing both hypersurfaces and illustrating their inter- 
section. Along this intersection curve of both hypersur- 
faces, \X\ (and so T*) is expected to be maximal for 
the CW signal examined in this example. The contour 
lines of constant / (and Z^^-*) in the sky of both hyper- 
surfaces are shown on the three-dimensional unit-sphere 
by Figure 3(d) and in a 2D Hammer- Aitoff projection by 



Figure 3(e) Here, in the sky subspace the intersection 
curve of Tii and 7^2 (corresponding to maximal detec- 
tion statistic) approximately coincides with the contours 
of Hz. 



D. Validity estimation of the global-correlation 
hypersurface description 

This section discusses how many hypersurfaces have to 
be considered in order to describe the detection-statistic 
maximum-structure. In other words, the question is 
investigated, in which region of the search parameter 
space and at which order m the power series (j29p can 
be truncated to approximate the matched-filtering am- 
plitude I A" I better than a certain fractional loss one tol- 
erates. 

For a given order m, one can estimate the contribu- 
tions from the next-order term in the matched-filtering 
amplitude \X\. As was done earlier in Section IIV CI 



here we consider again the four-dimensional parame- 
ter space {/, /^^^n} of "typical" wide-band all-sky CW 
searches, with / < 1 kHz, l/^^-*] < f/Tn^n, and a mini- 
mum spin-down age of Tmin = 50 yrs. 

To investigate the contribution from the first-order 
term in Equation ([50)1 to \X\, we compute the following 
integral ignoring terms in t with order m > 1, 



(42) 



Figure 4(a) shows |Xi | as a function of the dimensionless 



parameter Ami T. As already explained earlier, |Xi| at- 
tains its global maximum of one, when Aui = 0, which is 
the first order global-correlation equation. If one requires 
for instance |Xi| > 0.9, using Equation (|42)l this yields 
the condition |Ami|T < 1.57. When expressing Aui in 
terms of the original phase parameters using (fT9|) . one 
obtains 



lAuil < 2n 



|A/| 1+e 



|A/^'^IC + 2|/nC 



(1) 



(1) 



(43) 



But for typical wide-band all-sky CW searches with co- 
herent observation times T of a day or a few days, ac- 
cording to Equation |Ami|T can be considerably 
larger than 1.57. Therefore, it is clear that the first- 
order term in Equation (j30p contributes significantly to 
the matched-filtering amplitude \X\ and obviously can- 
not be neglected. 

Similarly, we estimate the contribution from the 
second-order term given that Aui = and by ignoring 
terms with m > 2 in Equation ([30]) . We calculate the 
corresponding integral denoted by X2 as 



X, . 1 
T 



T/2 



JAu2 r 



T/2 

e"/40F 

VAu2 



: erf 



dt 



(44) 



In Figure 4(b) IX2I is shown 0jS cL function of Au2 T . 
Obviously, IX2 attains the maximum of one at the sec- 
ond order global-correlation hypersurface, which is de- 
scribed by Au2 = 0. One finds that IX2I > 0.9 as long as 
IAM2I < 6.11. Reexpressing Am2 in the original phase 
parameters yields 



|Au2| < 27r 



(2) 



|A/ 



(1)1 



2 ^ 



(45) 

But according to Equation ([45]) for typical wideband 
all-sky CW searches with coherent observation times T 
of a day or a few days, |Au2|r can take values con- 
siderably larger than 6.11. Therefore, the second-order 
term in Equation (|30p also contributes significantly to the 
matched-filtering amplitude \X\ and cannot be neglected. 

Continuing this scheme to the next order in i, the con- 
tribution from the third-order term in Equation (|30p is 
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(a) Global-correlation hypcrsurfacc "Hi and contour lines of (b) Global-correlation hypcrsurfacc for 7^2 and contour lines 
constant frequency shown in the sky. of constant frequency shown in the sky (all close to each 

other). 




(c)Superposition of global-correlation hypersurfaces Hi^2i (d)Unit sphere with contour lines (circles) of 

and their frequency contour lines in the sky. global-correlation hypersurfaces Ti.i^2- 























-n/2 





RA Irad] 

(e) Hammer- Aitoff sky projection showing the contour lines 
(circles) of global-correlation hypersurfaces 7ii^2- 



FIG. 3: (Color online) The global-correlation hypersurfaces TCi (Awi = 0) and TC2 (Am2 = 0) for a given CW signal shown 



Hz/s. Each plot refers to GPS 



in the three-dimensional subspace {/, n} at the fixed target spin-down of j'^' = - 
time 793555944s. The phase parameters of the signal are as = 2.0 rad, 5s = — l.Orad 
indicated by the black cross (black circle) in the subspace {/, n} (in the sky plane). In (c) the superposition of (a 



-10" 

, /s 



100.0 Hz, /; 



(1) 



10-'"Hz/s, as 



and (b) 



illustrates the locations of expected maximum detection statistic along the intersection curve of both hypersurfaces. In each 
plot the light solid contour lines in the sky of constant / (and /'^') are of Hi, the dark solid contour lines are the ones 
of Ti.2- These contours of both hypersurfaces, shown in a Hammer- Aitoff projection in |(e)[ actually describe circles on the 
three-dimensional unit sphere as is visualized by (d) In the sky subspace of this example, the intersection curve approximately 



coincides with the contours of TI2, and due to the mismatch in spin-down the intersection curve does not pass through the 
signal's true sky-location. The black dotted-dashed curve represents the ecliptic. 
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considered. Given Atii = and Aw2 = 0, we estimate 
the contribution from the third-order term (with m = 3) 
by evaluating: 



1 

T 



T/2 

-T/2 
,^7^/6 



3 (AmsT^) 



3^l/3 



(46) 



where the gamma function r(a;) and the upper incom- 
plete gamma function r(a,x) are defined by 



r(x) 

r(a, x) 



i^-^e"* dt, 



-^-^e-Mt. 



(47) 
(48) 



As can also be seen from Figure 4(c) where Awa is 
plotted against jXaj, requiring jXal > 0.9 leads to the 
condition lAuaj < 9.79. By rewriting Au^ using again 
the original phase parameters according to (|19p , one ob- 
tains 



By approximating 

C(3) ^ o2 !^ 
S ~ "orb 



and 



(49) 



(50) 



a condition is obtained regarding the observation time T 
for which the contribution to \X\ from the third-order 
term in Equation (|30p is negligible: 



1/3 



(51) 



In the example presented in Section llV C[ a signal is con- 
sidered with /s = 100.0 Hz and /i^^ = 10-i°Hz/s. In this 
case, according to Equation ([STj) the third order can be 
neglected for observation times T < 61 h. For a similar 
signal and search, but for instance at /s = 1000.0 Hz the 
third-order term is expected to be non-significant only 
for observation times T < 28 h. 

However, considering cases where the third order in t is 
non-negligible, an analogous estimation of the contribu- 
tion from the fourth-order term is necessary. Therefore, 
we calculate X4 as 



X4 



T J -T/2 

2 



1/4 



(52) 
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FIG. 4: Estimating the number of global-correlation hyper- 
surfaces which contribute significantly to the detection statis- 
tic. In each plot both axes are dimensionless. Further details 
are given in the text. 



Figure [4(d)] shows 1X4 1 as a function of AU4T*. From 
the condition IX4I > 0.9 follows |Au4| < 27.78. Using 
this, analogously to Equations (|49p and ([^T]) . an approx- 
imate condition upon T can be found, which in the case 
of the present example signal yields that IX4I > 0.9 for T 
being approximately less than 10 days. As for the wide- 
band all-sky CW searches described earlier, coherent ob- 
servation times are typically only of one day or a few 
days, this means that in practice contributions from the 
fourth-order term in Equation (|30|) to \X\ are insignifi- 
cant in such cases. 



V. PREDICTIONS BY THE 
GLOBAL-CORRELATION EQUATIONS VERSUS 
FULL .7^-STATISTIC 



In this section the analytical predictions of the global- 
correlation equations ((35|) based on the simplified detec- 
tion statistic T* are compared to the results of fully co- 
herent searches using the detection statistic J- in different 
data sets containing artificial CW signals. The compu- 
tation of the full J-"-statistic includes effects of amplitude 
modulation and involves precise calculation of the detec- 
tor position with respect to the SSB using an accurate 
ephcmeris model. 
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TABLE I: Amplitude and phase parameters introduced 
in Section 1111 Bl of the two software-simulated continuous 
gravitational- wave signals. 



Signal 


A+ 




ip [rad] 


$0 [rad] 


1 


1.0 


0.5 


1.0 


2.0 


2 


1.0 


0.0 


0.0 


0.0 




Signal 


as [rad] 


6s [rad] 


h [Hz] 


fi''' [Hz/s] 


1 and 2 


2.0 


-1.0 


100.0 


-10-"' 



A. Comparison with the jT-statistic for 
software-simulated signals with no detector noise 

Different data sets have been prepared each contain- 
ing one of two different software-simulated CW signals 
without noise. A JT-statistic search has been conducted 
in each data set. The software tools used for data pro- 
duction as well as for data analysis are part of the LSC 
Algorithm Library Applications [l3| . 

The signal parameters defining the two simulated 
sources are given in Table HI In both cases the phase 
parameters are chosen to be identical. These are also 
the same phase parameters of the signal generating the 
hypcrsurfaces illustrated in Figure [31 Therefore, solely 
based on the global-correlation equations one expects the 
^-statistic to have a global large-value structure very 
similar to the one shown by Figure [3l To investigate 
the impact of the antenna pattern functions distinct am- 
plitude parameters have been chosen here, as given by 
TableHl For every set of simulated data the detector posi- 
tion refers to the LIGO Hanford 4-km (HI) detector, and 
the time of reference is chosen to be global positioning 
system GPS time 793555944 s consistent with Figure [3l 

As listed by Table HH for each data set containing 
the same signal different searches have been conducted, 
where the coherent observation time T has been var- 
ied between from 10 hours up to 2 sidereal days. In 
each search consisting of evaluating the JF-statistic on a 
grid of templates, an isotropic sky grid with equatorial 
spacing of 0.02 rad, spacings of 1/(2T) in the frequency- 
interval / € [99.8, 100.2] Hz and a fixed spin-down tem- 
plate of /'■^-' = -10^^^ Hz/s have been employed. 

First we compare the results of the searches with the 
prediction by the global-correlation Equation (|40ap form- 
ing the hypersurface Tii, shown in Figure 3(a) For each 
sky position and spin-down template-grid point where 
the JF-statistic search reported a candidate event with 
frequency /jf, one can calculate the relative deviation 
to the predicted frequency /jf* obtained from Equa- 
tion (|40ap by \fj^ — fy^* \l Jt* ■ For each search the maxi- 
mum relative deviations between predicted /f* and mea- 
sured fjr over the entire sky (and spin-down) are specified 
m TableUD and are of order 10 ^. The corresponding av- 
erage relative deviations are typically even one order of 
magnitude smaller. As in the simplified phase model (jl7p 



TABLE IL Comparison of predictions by the global- 
correlation equations based the simplified detection statis- 
tic J^* with the fully coherent .F-statistic search results 
using data sets containing software-simulated continuous 
gravitational-wave signals. The search labels, listed in the 
first column, containing the number 1 (number 2) refer to data 
sets containing only Signal 1 (Signal 2), whereas labels with 
different letters correspond to different observation times, as 
can be seen from the second column. From the obtained re- 
sults the maximum relative deviations from the predicted fre- 



quencies are specified in the third column. 


Search 


Coherent obser- 


Maximum deviation 


label 


vation time T 


1/^-/^*1//^* 


Al 


10 hours 


2.5 X 10'^ 


Bl 


1 sidereal day 


2.2 X 10"^ 


CI 


30 hours 


2.0 X 10"^ 


Dl 


2 sidereal days 


3.2 X 10"'^ 


A2 


10 hours 


2.6 X 10-*^ 


B2 


1 sidereal day 


2.4 X 10"^ 


C2 


30 hours 


2.3 X lO"*^ 


D2 


2 sidereal days 


3.5 X 10"'^ 



the Earth's spinning component has been neglected, here 
the observed frequency-deviations are consistent with the 
fact that the relative corrections to the frequency modu- 
lation originating from the Earth's spinning velocity are 
of magnitude Wspin/c ~ 10~^. 

Finally, FigureOpresents the results of each JT-statistic 
search projected on the sky. As we are only interested 
in the loudest J^-statistic values, only strong candidate 
events from the ones reported in the different searches 
A1,B1,C1 and Dl above a threshold of 2T > 2000 arc 
shown, and for the searches A2, B2, C2 and D2 a thresh- 
old of 2!F > 1250 is set. These thresholds were chosen 
to reduce the data volume to process and are based the 
expected global maximum values of 2J-'. These thresh- 
olds guarantee that candidate events with ^-statistic val- 
ues at least larger than w 15% of the global maximum 
are retained in each search. This thresholding is justi- 
fied because comparing candidate events with very small 
detection-statistic values to the global maximum are not 
of interest in this work. 

As shown earlier in Figure [31 in the two-dimensional 
sky projection the intersection curve of hypcrsurfaces Ti.i 
and Ti.2 approximately coincides with the contours (of 
constant / and fj^^^)of Ti2 in the sky. Thus, the dark 
circle in Figure 3(c) representing the contours of hy- 



persurface 7^2 also describes the predicted global max- 
imum structure of the detection statistic in the sky. In 
fact, this prediction is actually observed in qualitatively 
good agreement with the J-"-statistic search results shown 
in Figure [5] One finds that for coherent observation 
times less than one sidereal day (in searches Al and A2), 
the locations of the predicted global maximum struc- 
ture are only faintly visible in the results, because this 
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(a)Results of jF-statistic search Al (T = 10 h). 
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(c)Results of J'-'-statistic search Bl (T = Isd). 
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(e)Results of J'-'-statistic search CI (T = 30 h). 
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(g)Results of .7-"-statistic search Dl (T = 2sd). 
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(b)Rcsults of .7="-statistic search A2 (T = lOh). 
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(d)Rcsults of jF-statistic search B2 (T = Isd). 
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(f)Results of .F-statistic search C2 (T = 30 h). 
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(h)Resuhs of jF-statistic search D2 (T = 2sd). 



FIG. 5: Hammer- Aitoff sky projections of results from fully coherent jT-statistic searches in the data sets described by Table HIl 
Each data set contained one of the two software-simulated CW sources defined in Table H] where both signals 1 and 2 have 
identical phase parameters, but different amplitude parameters. The plots of the left (right) column show candidate events 
registered in the different searches A1,B1,C1 and Dl (A2,B2,C2 and D2) above a threshold of 2J" > 2000 {2T > 1250). The 
colorbar indicates the values of 2T x 10"'*. The corresponding analytical prediction by the global-correlation equations has 



been illustrated earlier by Figure 3(e 
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feature is still hidden due to the Earth's spin compo- 

For coherent observa- 



50 



nent [see Figures 5(a) and |(b)| 
tion times beyond one period of the Earth's spinning 
motion, the results clearly show the locations of large 

(h)]. As will be 



JF- values as predicted [see Figures 5(c 
discussed later in Section IVlIi the Earth's spinning mo- 
tion only varies the detection statistic within the global 
correlation-equations predicted (and observed) regions. 



B. Comparison with the jT-statistic for a 
hardware-injected signal in detector noise 

For so-called "hardware injections" simulated signals 
are physically added into the detector control systems 
to produce instrumental signals that are similar to those 
that are expected to be produced by astrophysical sources 
of gravitational waves. Through magnetic coil actuators 
the interferometer mirrors are made to physically move 
as if a gravitational wave was present. 

Here, we choose 30 hours of data (lying within a time 
span of less than 38 hours) containing a hardware injec- 
tion from LIGO's fourth science run (S4) of the LIGO 
Livingston 4-km (LI) detector. The GPS time of ref- 
erence is 795408715 s. During this data segment the 
hardware injection was activated 99.4% of the time. 
The signal's phase parameters arc defined by as = 
3.7579 rad, Ss ^ 0.0601 rad, /s 575.163636 Hz, and 



(1) _ 



— 1.37 X 10 ^^Hz/s, at the segment's starting 



time. The amplitude parameters are = 7.48 x 10~ , 
Ax = -7.46 X lO"^'^, ?A = -0.22 rad and the initial phase 
is $0 = 4.03 rad. In this segment of data an all-sky !F- 
statistic search has been carried out in the frequency in- 
terval of / G [575.048, 575.221] Hz, and in the spin-down 
range /(^^ G [-1.04 x 10-^1.04 x 10-i°]Hz/s. The 
grid of templates employed spacings of 5.70 x 10~® Hz 
in /-direction, separations of 3.23 x 10^^"Hz/s in J^^-*- 
direction, and an isotropic sky grid with equatorial spac- 
ing of 0.03 rad. 

In Figure [6] the results of the actual JF-statistic search 
are shown and compared to the prediction by the global- 
correlation hypersurfaces. As the search covers a range of 
spin-down templates there will distinct contours of hyper- 
surface 7^2 corresponding to each /"^^^-template. There- 
fore, the structure of maximum detection statistic J^* is 
expected to be an annulus in the sky referring to the 
magenta region in Figure 6(b) This annulus is framed 



by the different contour lines of Ti.2 which correspond to 
the minimum (dark-brown colored) and maximum (black 
colored) value of /^^^ searched. Qualitatively, the large- 
value structure of J- that is visible, is in good agree- 
ment with the predicted structure based on the global- 
correlation hypersurfaces. In Figure 6(a) the J^-statistic 
results only faintly reproduce the entire predicted struc- 
ture because of the Earth's spinning motion. Section IVIII 
will discuss and explain analytically how the Earth's 
spinning motion varies the detection statistic within the 
regions determined by the global-correlation hypersur- 
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(a)Rosults of a coherent all-sky J'-'-statistic search in a data 
set containing a hardware-injected signal. The colorbar 
indicates the values of 2JF. 
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(b)Prediction of the global maximum structure of the 
detection-statistic based on the global-correlation 
hypersurfaces. 



FIG. 6: Comparing the jT-statistic results (a) for a hardware- 
injected CW signal with the theoretical prediction (b) by the 
global-correlation equations. Both plots show Hammer-Aitoff 
projections of the sky. The sky location of the simulated 



signal is represented by the black cross in (b) The magenta 
region in |(b)| represents the predicted structure of maximum 
detection statistic J^* and is observed to agree qualitatively 



well with the actual jT-statistic results (a) from the hardware 
injection in real detector data. 



faces. 



VI. VETOING INSTRUMENTAL NOISE 
ARTIFACTS 

A typical feature of the data from an interferometric 
gravitational-wave detector are narrow-band noise arti- 
facts, so-called "lines" , which are of instrumental origin. 
As a consequence, the results of a search for continuous 
gravitational-wave sources contain instrumental artifacts 
that in some respects mimic CW signals. But these arti- 
facts tend to cluster in certain regions of parameter space. 
For the case of incoherent searches as reported in (l3| . 
candidate events in such parameter-space regions were 
identified and discriminated. Here, we propose a simi- 
lar veto method also suitable for coherent CW searches, 
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such as those found in Using the global-correlation 
equations found in this work we aim to identify those re- 
gions in parameter space where instrumental noise lines 
can imitate a real signal by producing large detection- 
statistic values. In such a case these candidate events 
could automatically be vetoed. 

For simplicity we consider the same four-dimensional 
parameter space as used in Section IIVBI consisting of 
{/, f^^\ n} and that T has a value, such that third-order 
contributions to \X\ arc insignificant (cf. Section riVD|) . 
Thus, the global-correlation equations of relevance are 
the same as in the example studied in Section IIV CI 
These were given by Equations (|40a[) and (j40bp . It is 
obvious that the frequency of a stationary instrumental 
line is independent of the Earth's position in its orbit 
around the sun. This decoupling is achieved by setting 
Hs = 0. For a stationary instrumental line originating 
from the detector it also holds that /s^-* = 0. In this case 
the constants Ki^2 in Equations (|40ap and (|40bp simplify 
to Ki = /s and K2 = 0. Thus, the two relevant global- 
correlation equations are of the following form, 



;(i)+/^(2) 



n 



(53a) 
(53b) 



As stated earlier, for a given f'-'^^ Equation (|53ap de- 
scribes a three-dimensional hypersurface in the subspace 
{/, n}. On this hypersurface the detection statistic 
will attain its maximum along the intersection curve 
with the hypersurface described by Equation (|53bp . In 
a projection into the sky subspace this intersection 
curve approximately coincides with the contours (of con- 
stant / and Z^^-*) of hypersurface (j53bp . Therefore Equa- 
tion (j53bp describes the region in the sky for which po- 
tential CW signals do not produce a modulation pattern 
that would distinguish them from an instrumental line. 

Using this knowledge one can discriminate (veto) can- 
didate events which satisfy Equation (j53bp . As the reso- 
lution in parameter space is finite, we formulate the fol- 
lowing veto condition: 



(2) 



n 



(1) 



n 



(54) 



the tolerance-parameter e > can be understood as 



AT 



coll; 



(55) 



where A/ccii denotes the resolution in the frecj^ucncy- 
direction (width of cells), A^coii the number of cells one 
tolerates during a characteristic length of time AT. 

One can visualize and calculate the volume of the re- 
gion in four-dimensional parameter space which is ex- 
cluded by this veto. For a given source sky position. 
Equation ([5^ is linear in / and f^^\ Thus, for fixed sky 
position n, the veto condition defines two parallel lines in 
the {/, Z^^-*} — plane. Candidate events which lie in the 
region between the lines are discarded (vetoed). Can- 
didate events which lie outside this region are retained 



(not vetoed). The locations of these two lines in the 
{/j/^^^} — plane depends upon the sky position. The 
fractional volume excluded by the veto depends upon 
whether or not (as the source position varies over the 
sky) the excluded region between the lines lies inside or 
outside of the boimdaries of the search, or intersects it. 
Alternatively, for a given value of / and f'-'^\ one can 
calculate the portion of the sky which is excluded by 
the veto, depending upon the ranges of parameter space 
searched. The details of such a calculation for a particu- 
lar search can be found in Appendix A of . 

Figure [7| illustrates the veto method for an example 
noise line. Thereby, Figure 7(a) presents the results of 
a fully coherent matched-filtering search using the JF- 
statistic for a 30-hour observation time. The data set an- 
alyzed contains a detector-noise line, which is in this case 
a violin mode resonance of the mode cleaner mirrors of 
the LIGO Hanford 4-km (HI) detector. In Figure |7(b)"| a 
comparison with the theoretical prediction by the global- 
correlation equations given by (|53ap and (|53b[) is made 
featuring a very good agreement. Thus the veto will be 
very efficient in excluding the parameter-space regions of 
largest 2J^-values produced by instrumental lines. 

Note that this veto excludes only the loudest candi- 
date events (of largest 2!F) , but as shown in Figure [7] an 
instrumental line is capable of contaminating large parts 
of the sky due to the global correlations (depending on 
the search parameters) . Therefore in some cases it might 
be necessary to increase the tolerance-parameter e artifi- 
cially to account for these effects if necessary. However, in 
a J^-statistic search one is interested in the strongest can- 
didate events arising from the background level. Thus in 
eliminating such false instrumental-noise events the veto 
condition presented is very efficient as these regions are 
well described. This veto method as presented here is 
applied in the EinsteinQHome search described in [13] . 



VII. EFFECTS OF THE DIURNAL SPINNING 
MOTION OF THE EARTH 



By considering the Earth's spinning motion in the 
phase model, we here investigate the variation of the 
detection statistic along the predicted global maximum 
structure by the global-correlation hypersurfaces. In 
other words, given the solution Aw,„ = 0, which are the 
the global-correlation hypersurfaces, we study how the 
spin component 0spin(i) of Equation (jl3p in the phase 
model ([4]) modulates the detection statistic in the lo- 
cations consistent with this solution. In order to sim- 
plify this discussion, in what follows frequency time- 
derivatives in the spin component </'spin(0 ignored 
to obtain 



nit) ■ n 



(56) 
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where we defined the vector 
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(a)Rcsults of a coherent all-sky jF-statistic search in a data 
set containing an instrumental noise line. The colorbar 
indicates the values of 2T. 
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(b) Prediction of the global maximum structure of the 
detection statistic based on the global-correlation 
hypersurfaces. 

FIG. 7: Comparing the results of a fully coherent matched- 
filtering search using the .7^-statistic (a) with the theoretical 
prediction (b) by the global-correlation equations for a given 
instrumental-noise feature mimicking a real CW signal. Both 
plots show Hammer-Aitoff projections of the sky. The all- 
sky search was carried out in a 0.5 Hz frequency-band / G 
[568.0, 568.5] Hz and for a range of frequency time-derivatives 
/'^^ G [-3.63 X 10~^3.63 x 10"^"]Hz/s, for an observation 
time of T = 30 h. The GPS time of reference is 795149016 s. 



The upper plot (a) shows all candidate events reported by the 
search above a detection-statistic threshold of approximately 
50% of the largest 2.7^-value found. The frequency of the 
instrumental-noise line present in this data set is a resonance 
violin mode of the mode cleaner mirrors of the LIGO Hanford 
4-km (HI) detector. The magenta region in |(b)| of maximum 
expected detection statistic agrees well with the .7^-statistic 
results (a) from the real instrumental line. 



Thus, the phase difference between the spin component 
of the signal (/)spjii(i) and a template 0spin(i) is given by 



Ak = 27r (ns/s - n/) 



(58) 



Taking into account A(/)spin(t) in the detection-statistic 
amplitude and provided that Au„i = 0, one has to com- 
pute the following integral 



1 

T 



T/2 



-T/2 



(59) 



For observation times T relevant to CW searches (of 
order days), the phase modulation due to the spinning 
motion of the Earth is oscillatory, because it has a period 
of one sidereal day figpin = 27r/lsd. Therefore, in order 
to evaluate ((5^ we follow a route previously taken in 0, 
[isl , which makes use of the Jacobi- Anger identity: 



E 

n— — oo 



^"J„(2)e'' 



(60) 



where Jn{z) is the n-tli Bessel function of the first kind. 
This identity allows to expand exponentials of trigono- 
metric functions in the basis of their harmonics. To 
employ the Jacobi-Anger expansion we rewrite Equa- 
tion (|57p by approximating the diurnal detector motion 
due to the Earth's rotation to be circular. 



A0spi„(<) = 



A/cii sin A 



-|-Afc^ cos A COs{(f)Q -f f^spini) 



(61) 



where i?E is the radius of the Earth, A is the latitude of 
the detector, Afcy is the absolute value of the component 
of the vector Ak parallel to the rotation axis, A/c^ is 
the absolute value of the component of Ak orthogonal to 
the rotation axis, and (/jq is determined by Ak at t = 0. 
Defining 



-Re 



spin, 



A/cii sin A, 



(62) 
(63) 



Augpin = — A/cj^ COS A, 
c 

we apply the Jacobi- Anger identity: 

oo 

g»A03p.„(t) ^j,»A03p.„,ii i"e"'^''J„(Auspi„)c'' 

n— — oo 

(64) 

Substituting this expression into Equation (|59p and tak- 
ing the modulus, one obtains 



S D i 11 ^ 



A0spin(t) 



0spin(O - 0spin(O 
I'spin (^) 



27r 



c 

I'spin (^) 



(ns/s - n/) 



Ak, 



(57) 



l-'i^spinl 



f; z"c-^V„(Au,pi„)sincf^%^ 



n— — OO 



(65) 

Figure [5] shows jXspinj for the two LIGO detectors as 
a function of T and Augpin. It is obvious that for the ob- 
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(a)LIGO Hanford detector 
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(b)LIGO Livingston detector 

FIG. 8: (Color online) The simplified detection-statistic am- 
plitude |Xspin| for phase- mismatch only in the spin compo- 
nent of the phase model as a function of observation time T 
and dimensionless parameter-mismatch Aitspin as defined in 



Equation (|63|l . for the two LIGO detectors. For observation 
times beyond T > 27r/r2spin, a good approximation of |Xspin 
is given by the dominant term | Jo(Auspin)|. 



servation time being an integer multiple £ of the Earth's 
spin period, such that T = 27r£/57spin, Equation fSS]) 
then simplifies to |Xspin| = | Jo(Auspin)|, because only 
the term corresponding to n — does not vanish. By in- 
spection, we find that this is also approximately the case 
for all observation times of one day or longer, as can be 
seen from Figure [H] Therefore, we approximate |Xspin| 
foi' ^ 27r/rispin (which is also the regime relevant to 
CW searches) by 



\X, 



spin I 



I Jo(AWspin)|- 



(66) 



Figure [ni illustrates |Jo(Auspin)| over the entire sky for 
the three different cases studied previously in this work: 
the software-injected signal, the hardware-injected sig- 
nal and the instrumental-noise line. Comparing Fig- 
9(a) with Figure [5l Figure 9(b) with Figure 6(a)[ 



ure 



and Figure 9(c) with Figure 7(a)[ one finds that the 



variation of the J^-statistic in the regions determined by 
the global-correlation hypersurfaccs (locations consistent 
with Au„i = 0) can in fact be recovered in Figure [5) 
To illustrate better why the Earth's orbital motion de- 



termines the regions of largest detection-statistic and the 
Earth's spinning motion only modulates the detection- 
statistic within these regions, we consider the examples 
represented in Figure [TUl There, as a basis of comparison 
to |Xspin| of Equation ((66)l we use \X\ of Equation pS]) 
of the orbital phase model for m < 2. As shown in Sec- 
tion IIVDI higher orders of m are insignificant for the 
example cases considered. Figure [TU] compares \X\ and 
l-'^^spinl as functions of sky position for the given signal 



phase parameters. In the two diagrams 10(a) and (c) 



the sky has been sliced along a at constant declination 
(5 = (5s, and for fixed values of frequency / = /s and for 



simplicity also at zero spin-down offset /< 



(1) 



In 



the plots |10(b)| and |(d)[ the sky has been sliced along 
declination 5 at constant right ascension a = as-, and for 
the remaining template parameters coinciding with the 
signal's parameters. The essential observation is that 
\X\ drops off much more rapidly compared to |Xspi,i|. 
Therefore, \X\ dominantly determines the global maxi- 
mum structure of the detection statistic, whereas |Xspin| 
only modulates the detection statistic within these re- 
gions where \X\ is maximal. 



VIII. CONCLUSIONS 

The family of global-correlation hypersurfaccs derived 
in this paper provides a approximate analytical descrip- 
tion of the global large- value structure of the optimal de- 
tection statistic J-^ in the phase-parameter space of con- 
tinuous gravitational-wave searches. 

For observation times longer than one sidereal day, but 
still much smaller compared to one year, it is the orbital 
motion of the Earth which generates a family of global- 
correlation equations. The solution to each of these equa- 
tions is a different hypersurface in parameter space. The 
detection statistic is expected to have large values at the 
intersection of these hypersurfaccs. In this context, the 
Earth's spinning motion plays a minor role, because it 
only varies the detection statistic within the intersection 
regions determined by the global-correlation hypersur- 
faccs. 

While embedding previously published results Q in 
the present theory, this work leads to a substantially im- 
proved understanding of the global correlations in the 
optimal detection statistic. 

In a comparison study with results of the J-'-statistic 
from numerically simulated as well as from hardware- 
simulated signals in the presence of noise, the analytical 
predictions by the global-correlation equations have been 
qualitatively well recovered. 

In addition, the global large- value structure of the de- 
tection statistic produced by stationary instrumental- 
noise lines mimicking astrophysical sources is also well 
described by the global-correlation equations. This per- 
mitted the construction of a veto method, where such 
false candidate events are excluded. 

Moreover, reparameterization of the original phase 
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FIG. 9: (Color online) Shown is |Xspin| ~ | Jo(Auspin)| over the entire sky using equatorial coordinates of right ascension a and 
declination 5 for the three different examples studied previously in this work. The first column corresponds to the software- 
simulated signal [cf. Figure [5] 
the instrumental-noise line [cf. Figure 7(a 



the second column to the hardware-injected signal 
In the top row. 



cf. Figure 6(a 
three-dimensional plots of 



and the third column to 
as functions of a and 5 



Jo(AUspm 

are shown, and in the bottom row two-dimensional Hammer- Aitoff projections of the sky are given illustrating the contours of 
I Jo(Atispi 



parameters by the parameters Um from the global- 
correlation equations offers evident advantages in solving 
further problems related to CW searches. For example 
using the Wm-parameters can help in placing templates 
more efficiently. As these parameters "absorb" the global 
correlations leading to a linear phase model, the metric in 
these parameters will be flat, which means independent 
of the parameters (cf. [3). 

Furthermore, one could implement a candidate-event 
coincidence scheme via reparamcterizing the phase pa- 
rameters of the candidate events by the Um-parameters. 
In the results from a coherent search using a bank of 
templates, a putative signal should accumulate large de- 
tection statistic values in a very focused region in the 
space of Mm-parameters, because for the templates lo- 
cated along the contours of the global large- value struc- 
ture, the parameters Um are approximately invariant. In 
addition this fact gives rise to a novel type of hierar- 
chical semi-coherent search technique for CW sources. 
In such a multistage scheme d, [l3| one breaks up the 
data set into a sequence of short data segments, of which 
each segment is analyzed coherently in a first stage. This 



is followed by an incoherent combination of the coher- 
ent results from each segment. The transformation to 
global-correlation parameters Um helps classifying coin- 
cident candidate events from the first-stage coherent step 
before these are combined incoherently. However, in fu- 
ture work the efficiency of such a correlation transform 
scheme [Ol as well as the relevance of the global corre- 
lations in the context of the template-placing problem 
should be investigated. 
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l^apinl are shown as functions of declination template- value 5, while the remaining template parameters coincide with the 
signal's phase parameters. One can see that \X\ decreases much more rapidly compared to |Xspin|. Therefore \X\ dominantly 
determines the global maximum structure of the detection statistic. 
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